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ABSTRACT 



Our recent progress on time-dependent modeling 
of the multiwavelength spectra and variability of 
blazars with leptonic and hadronic jet models is re- 
viewed. Special emphasis is placed on X-ray spectral 
variability of low-frequency peaked (LBLs) and in- 
termediate BL Lac objects (IBLs). As an example, 
recent observational and modeling results of an ex- 
tensive multiwavelength campaign of BL Lacertae in 
2000 are presented. It is demonstrated how com- 
bined spectral and variability modeling of LBLs and 
IBLs can significantly constrain emission models and 
potential variability scenarios. In the case of BL Lac- 
ertae, the variability appears to be driven primar- 
ily by fluctuations in the spectral index of the non- 
thermal, ultrarelativistic electron population in the 
jet. Such constraints allow us to refine predictions of 
the intrinsic GeV 7-ray emission and the dominant 
electron cooling mechanism in these objects. 

Key words: Active galactic nuclei; blazars; BL Lac- 
ertae; jets; gamma-rays; multiwavelength observa- 
tions; theory. 



1. INTRODUCTION 



Blazars are the most extreme class of active galaxies 
known to date. They have been observed in all 
wavelength bands — from radio through very-high 
energy (VHE) 7-ray frequencies. 66 blazars have 
been identified as sources of > 100 MeV emis- 
sion detected by the EGRET telescope on board 
the Comvton Gamma-Ray Observatory (CGR'^ 
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at VHE 7-rays (> 350 GeV) by ground-based air 
Cerenkov telescopes. Many of the EGRET-detected 



7-ray blazars appear to emit — at least temporar- 
ily — the bulk of their bolometric luminosity at 
7-ray energies. Blazars exhibit variability at all 
wavelengths on time scales — in some cases — 
down to less than an hour. VLBI radio observations 
and monitoring often reveal one-sided kpc-scale 
jet structures, exhibiting apparent super luminal 
motion. The radio through optical emission from 
blazars often shows linear polarization, pointing 
towards a synchrotron origin. 



1.1. Blazar Spectra 

The broadband continuum spectra of blazars are 
dominated by non-thermal emission and consist of at 
least two clearly distinct, broad spectral components. 
A sequence of sub-classes of blazars can be defined 
through the peak frequencies and relative vF^ peak 
fiuxes of those components, which also appear to be 
correlated with the overall bolo metric luminosity of 
the sources l|Fossati et al.lll998D : In the case of flat- 
spectrum radio quasars (FSRQs), the low- frequency 
(synchrotron) component extends from radio to opti- 
cal/UV frequencies, with a peak frequency generally 
in the mm or IR band; the high-frequency compo- 
nent extends from X-rays through GeV 7-ray ener- 
gies, with a vFn peak frequency corresponding to 
~ 10 MeV — 1 GeV. No FSRQ has so far been de- 
tected by ground-based air Cerenkov telescope fa- 
cilities at energies > 100 GeV, although in flaring 
states the 7-ray vF,j peak flux of FSRQs dominates 
over the low-frequency emission by up to ~ 1 order 
of magnitude. In the case of high-frequency peaked 
BL Lac objects (HBLs), the low- frequency compo- 
nent often extends far into the X-rays, with peak 
frequencies ranging from the U V/soft X-ray to the 
hard X-ray regime ijPian et al.l[l998.) . depending on 
the source and its state of activity; the high-energy 
component of HBLs extends from hard X-rays far 
into the VHE 7-ray regime. All blazars detected 
at VHE 7-ray energies to date are HBLs. In spite 
of extending to extremely high photon energies, the 
vF^ peak flux of the 7-ray component of HBLs is 
generally at most comparable to the spectral output 
in the low- frequency component. In terms of their 
overall bolometric luminosity, FSRQs appear to be 
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Figure 1. Broadband spectra of th e FSRQ 3C 279 (left y anel: \Hartma,n et a,l\ . \200h\:\CoU,m.ar et al . \200A ). the 
LBL BL Lacertae (middle vanel ABottcher et ai\.\200^) . and the HBL Mrk 501 f Petrv et al\. \20m} . For each 
object, two simultaneous SEDs at two different epochs are shown, illustrating the range of long-term spectral 
variability see n in these objects. Fo r 3C 279, the results of recent multiwavelength observations in 2003, including 
INTEGRAL \Collm,a.r et alV^EoI) are also shown. 



several orders of magnitude more powerful sources 
than HBLs. 

Apparently intermediate between the FSRQs and the 
HBLs are the low-frequency peaked BL Lac objects 
(LBLs). The peak of their low- frequency component 
is typically located at IR or optical wavelengths, their 
high-frequency component peaks typically at ~ sev- 
eral GeV, and the 7-ray output is of the order of 
or slightly higher than the level of the low-frequency 
emission. Fig. ^ shows a compilation of broadband 
spectra of a typical FSRQ (3C279), an LBL (BL Lac- 
ertae), and an HBL (Mrk 501), respectively. Each 
panel shows the respective object at two or three 
different observing epochs, in two different activity 
states. 



1.2. Blazar Variability 

Fig. n already illustrates that in particular the high- 
energy emission from blazars can easily vary by more 
than an order of magnitude between different observ- 
ing epochs. However, variability has been observed 
on much shorter t ime scales, in some extreme cases 
less than an hour ijGaidos et a l.. 1996). Fig. |21shows 
examples of light curves of the LBL BL Lacertae, 
taken during a broadband observing cam paign car- 
ried out in the seco nd half of 2000 (Bott cher et al.l 
l2005lVillata et al.ir20021 . While the radio emission 
of blazars is generally variable on time scales of weeks 
- months (Fig. |2K), the optical light curve shows sig- 
nificant variability on time scales down to ~ 1.5 hr 
(Fig. Eb). 

Often, both the optical and X-ray emission show 
characteristic hardness- intensity correlations. Fig. 
O illustrates this for BL Lacertae in 2000. Some 
HBLs (e.g., Mrk 421 and PKS 2155-304) have been 
observed to exhibit characteristic , clockwise loop 
struc t ures ( "spectral hyster esis" ; iTakahashi et al.l 
IT99I iKataoka et all l2000ll . which can be in- 
terpreted as the synchrotron radiation signature 
of gradual injection and/or acceleration of ul- 



trarelativistic electrons into the emitting region, 
and subsequent radiative coo l ing feirk et all 119981: 

]1998f; 'Kataoka et aL|, 
Li fc Kusunosc, 2000). 



|Georganopoulos_^^fe|rscliei 
iboO; Kusunos Tet alJ.l2000t 



In LBLs, the soft X-ray emission is also sometimes 
dominated by the high-energy end of the synchrotron 
emission component, so similar spectral hysteresis 
phenomena should in principle be observable. How- 
ever, those objects are generally much fainter at X- 
ray energies than their high-frequency peaked coun- 
terparts, making the extraction of time-dependent 
spectral information an observationally very chal- 
lenging task. Fig. 121 clearly illustrates that the 
BeppoSAX observations of BL Lacertae in 2000 re- 
vealed evidence for spectral variability, but lacked 
the sensitivity to clearly establish or rule out spec- 
tral hysteresis. Such a measurement might require 
the new generation of X-ray telescopes such as Chan- 
dra or XMM-Newton. 



OVERVIEW OF LEPTONIC JET MODELS 
OF BLAZARS 



The high apparent bolometric luminosity combined 
with the short variability time scales and the appar- 
ent superluminal motions of individual jet compo- 
nents observed in many blazars, provide compelling 
evidence that the nonthermal continuum is produced 
in emission regions of a typical size scale of ~ a few 
light days or less, moving relativistically along a jet 
structure which is directed at a small angle with re- 
spect to our line of sight. The jets are most likely 
powered by accretion of circumnuclear matter onto a 
supermassive black hole of 10^ Mq <Mbh < 10^° Mq. 
The emission regions are characterized by the pres- 
ence of an ultrarelativistic population of nonthermal 
electrons. Several scenarios have been proposed con- 
cerning the acceleration of such ultrarelativistic elec- 
trons, including impulsive injection near the base of 
the je t Ipermcr & Schlickciser, 1993; Bottcher et aQ, 
Il997j) . individual shock waves propagating along 
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Figure 2. Radio and optical light curves of the LBL BL Lacertae during the multiwavelength campaign of 2000. 
The dark- shaded areas in panel a) indicate the time segments of BeppoSAX observations. 



the jet iMarscher fc Gean l)1985fl . relativistic par- 
ticle acceleration at shear layers between a fast- 
moving inner jet and a slow er moving outer jet (e.g., 
iStawarz fc Ostrowski 12003(1 or internal shocks from 
the collisions of multiple shells of material ejected 
into the jet structure (SDadactal., 200^. Be- 
cause of the difficulty of constraining the acceleration 
mechanism and the composition and spectral char- 
acte ristics of the injected particle distribution (see , 
e^..'Sikora fc Madeiski'.'200d':'Ost rowski fc Bednard. 
12002: Stawarz fc Ostrowski, 2003J), the time profile 
of injection and injected particle spectra of ultrarel- 
ativistic electrons are generally treated as free pa- 
rameters in blazar modeling. 

The nonthermal electrons are emitting synchrotron 
radiation, which is responsible for the low-frequency 
emission from radio to UV or even X-ray frequen- 
cies. Higher-frequency (X-ray and 7-ray) emis- 
sion is produced via Compton scattering processes. 
Possible target photon fields for Compton scatter- 
ing are the synchrotron p hotons produced within 
the jet (the SS C proce ss: IMarscher fc Geail 1 19851 
iMaraschi etldl 1992: Bl oom fc MarscheA 11996(1 . or 
external photons (the EC process) . Sources of exter- 
nal seed photons include the UV - soft X-ray emis- 
sion from the disk — - either entering the jet directly 
l|Dermer et all Il992t iDermer fc Schlickeiseil Il993(l 
or after reprocessing in the broa d line region ( BLR) 
or other circumnuclear material ('Siko ra et ali ll994t 
iBlandford fc Levinsod ri995,: Dermer et all Il997|) 
— , j et synchrotron rad i a.tion reflected a t the 
BLR dC hisellini fc Madau', 1996'; 'Bednarek", 1998^ 
iBottcher fc Dermer, 1998), or the infrared emission 
from circumnu clear dust (B lazeiowski et al.L 120001: 
lArbeiter et all [2002) . In the context of jet model 
invoking a significant deceleration of the outflow 
along the jet, an important source of soft photons 
might also be provided by the synchrotron photons 
from slower portions of the jet further downstream 



l|Georganopoulos fc Kazanaa 12003(1 . 

In addition to these fundamental radiation processes, 
77 absorption and pair production as well as syn- 
chrotron self absorption have to be taken into ac- 
count in order to build a self-consistent blazar radia- 
tion model. Synchrotron self absorption is the reason 
why some of the spectral fits shown in Fig. ^do not 
reproduce the radio spectra of the respective sources: 
The emission regions are optically thick at radio fre- 
quencies during the early stages of propagation of 
the emission regions through the jet, where the non- 
thermal electron population is sufficiently energetic 
to produce ample high-energy radiation. 

As the emission region is moving relativistically out- 
ward, the nonthermal electron population will evolve 
according to 



dt 




ne[l,t] 



acc/loss,cont. 



+ 9e(7, t) - Veil, t) + Qe(7, 



nc{e,t) 



(1) 



where (d7/(ii)acc/ioss,cont. denotes the continuous en- 
ergy losses due to radiative and adiabatic cooling 
and energy gain due to acceleration processes, pe and 
qe are the terms describing the population and de- 
population of a given electron energy interval due to 
non-continuous energy loss processes (such as Comp- 
ton scattering, in particular in the Klein-Nishina 
regime), Qe describes the electron injection function, 
and io.csc is the escape time scale of nonthermal elec- 
trons. The evolution of the photon population in the 
emission region has to be solved simultaneously with 
the electron distribution, and is determined through 
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Figure 3. Left panel (a): Optical color-magnitude relation observed during the multiwavelength campaign on BL 
Lacertae in 2000 (Villa ta et al., 200k; Bott cher et al, , 2005; B ditcher &^ R,eim,er , 2001). Right panel (b): X-ray 
hardness-intensity diagram extracted from BeppoSAX observat i ons. T he curves indicate the spectral variability 
patterns from the best- fit leptonic model of \Bdttcher & ReimeA \200A) . 



dnpi,{e,t) 
dt 



^ph,cm ^ph,abs 

nphie,t) 



+ gph(e,t) -_pph(e,t) 



(2) 



where now nph,Gm and riph.abs describe the funda- 
mental emission and absorption processes, Qph and 
Pph describe the scattering rates into and out of a 
given photon energy bin, and tph.csc is the photon 
escape time scale. 



2.1. Leptonic- Jet Spectral Modeling Results for 
Different Blazar Classes 



Various versions of the generic leptonic jet model 
described in the previous section have been used 
very successfully to model simultaneous broadband 
spectra of several FSRQs, LBLs, and HBLs. As 
more detailed spectral information has become avail- 
able, the results of such broadband spectral model- 
ing have now converged towards a rather consistent 
picture l|Ghisellini et all 119981: iKiibo et all ^9§). 
The spectral sequence HBLs LBLs FSRQs 
appears to be related to an increasing contribution 
of the external Comptonization mechanism to the 
7-ray spectrum. While most FSRQs are success- 
fully mo delled with external Comptonization mod- 



els ('e.g..lD^mer et a].Ln997HSa,mbnina et a,lJ.I1997t 
iMukheriee et all [1999; Hart man et all l2001afl . the 



broadband spectra of HBLs are consistent wit h pure 
SSC models ( Mastichiadis & Kirk, 1997; Pia n et all 
liM IPetrv et a,U QMi IKrawczvnski et all 12000). 
BL Lacertae, a LBL, appears to be intermediate 
between these two extremes, requiring an external 
Comptonization component to explain the EGRET 
s pect rum (Madeiski et al., 1999; Bottcher & Bloom, 
l200rih . One generally finds that HBLs require higher 



average electron energies and lower magnetic fields 
than LBLs and FSRQs. In most cases, the required 
Doppler boosting factors D seem to be comparable 
for all types of objects, although there have also been 
some results indicating an extraordinar ily high value 
of_D>40 for the HBL Mrk 501 (Krawczvnski et all 
|2002() . Typical examples of broadband spectral fits 
consistent with this sequence of leptonic jet model 
parameters are shown in Fig. ^ The occasional find- 
ing of very high Doppler factors, in particular in some 
HBLs, has prompted iGeorganopoulos fc Kazana^ 
((2003) to propose their model of a decelerating, strat- 
ified jet in which synchrotron photons from slower 
regions of the jet would serve as seed photons for 
Compton scattering, appearing slightly blue shifted 
in the rest frame of the faster high-energy emission 
region further upstream. Such a scenario could re- 
move the need for bulk Lorentz factors largely in ex- 
cess of ~ 10 for HBLs. 



2.2. Leptonic- Jet Modeling of Blazar Spectral 
Variability 

The generic blazar model described above is inher- 
ently time-dependent and facilitates the modeling 
not only of the broadband SEDs, but also the 
detailed spectral variability of blazars. Studies of 
blazar spectral variability have been done in great 
detail for the case of pure SSC models with electron 
cooling dq mii iated by synchrotron losses, e.g., by 
iKirketal] ||1998D: IGeorganopoul os fc M arschei 
l|1998[) : iKataoka et al.l ll2000D : [kusunos e et al. 
l|2000() : iLi fc Kusunosd 1)20001) . In those papers, 
the spectral hysteresis observed in several HBLs 
was reproduced, significantly constraining model 
parameters beyond constraints obta inable from 
pure s pectral modeling. More recently, Sikor a et al.l 
l)2001() have extended these studies to an inho- 
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Figure 4- Comparison of the tracks in the flux - 
spectral-index plane for different generic jet models, 
from synchrotron- dominated {sh ort- dashed) to SSC 
domin ated (long- dashed). From iBottcher & Chiang 
MM) . 



iBottcher fc Chiand lj2002t) . A key result of their 
study was that a dominant SSC component would 
leave very obvious imprints in the X-ray spectral 
variability of these objects, as illustrated in Fig. 
^ In contrast, they found that a moderate and 
even slightly dominant contribution from external 
Compton scattering will have virtually no effect 
on the X-ray spectral variability patterns. This 
might be a consequence of the fact that the beaming 
pattern of external Commpton emission is more 
strongly peaked in the forward direction than 
the synchrotron and SSC components (which are 
assumed to be isotropic in the co-moving frame 
of the emission region). Thus, even if the EC 
component is dominating the 7-ray emission at 
MeV - GeV energies, it may only make a moderate 
co ntribution to the electron c ooling rate. The results 
of IBottcher fc Chiand ^^2002^ have been apphed in 
detail to simultaneous m ultiwavelength o bserv ations 
of W Comae in 1998 (Bottchcr et al., 2002) and 
BL Lacertae in 2000 (,Bottcher fc Rcimer. .2004) . 



mogeneous jet model, also including an external 
Compton component, so that the model would now 
also be applicable to LBLs and quasars. They have 
applied their model to the special case of 3C 279. 
Based on their results, they interpreted the lack of 
a measurable time lag be tween the 7-ray and X- ray 
emission in that object l|Hartman "etail l2001bl) as 
evidence that X-rays and 7-rays might be produced 
co-spatially by electrons of similar energies. This, 
in turn, provides evidence for two separate emission 
components being dominant at X-rays and 7-rays. 
Most plausibly, this might indicate that the X-ray 
emission is dominated by SSC emission, while 
the 7-rays are dominated by external Compton 
emission. 

Time-dependent, homogeneous leptonic jet models in 
which radiation mechanisms other than synchrotron 
may dominate have recently been investigated 
in an analytical approach by IChiang fc Bottcheil 
ll2002D. and with det ailed nu merical simulations 
bv 'Bottche r fc Chiand l|2002D . IChiang fc Bottched 
(|2002) pointed out that a dominant contribution 
from SSC to the electron cooling will produce a char- 
acteristic time-averaged synchrotron spectral index 
of a ~ 3/2 (energy spectral index), independent of 
the injection index of relativistic electrons in the jet. 
In this case, the cooling time scale of electrons radi- 
ating at a synchrotron photon energy Esy is expected 

to scale as Tcoo1,ssc oc Ei§ ^^1"^ ^ where q is the elec- 
tron injection spectral index (Bottchcr et al., 2003). 
This differs characteristically from the synchrotron 

— 1/2 

or EC dominated case in which Tcooi.sy cx; E^y , in- 
dependent of the electron injection index q. 

Detailed numerical simulations of the time- 
dependent emission characteristics of homogeneous 
jet models with parameters specifically chosen 
to be appropriate for low-frequency peaked and 
intermediate BL Lac objects have been done by 



3. OVERVIEW OF HADRONIC BLAZAR 
MODELS 



While leptonic models deal with a relativistic e 
plasma in the jet, in hadronic models the relativis- 
tic jet consists of a relativistic proton (p) and elec- 
tron (e~) component. In the following, a brief 
summar y of the Synchrotro n-Proton Blazar (SPB- 
) model l|Miicke et ali l2003^ is given, as an example 
of a hadronic model that takes into account all the 
salient features of hadronic blazar jet models in gen- 
eral. 

Like in the leptonic model, the emission region in 
an AGN jet moves relativistically along the jet axis 
which is closely aligned with our line of sight. Rela- 
tivistic protons, whose particle density rip follows a 
power law spectrum oc 7p '^'^ in the range 2 < 7^ < 
7p,max, are injected instantaneously into a highly 
magnetized environment (B = const, within the 
emission region), and are subject to energy losses 
due to proton-photon interactions (meson produc- 
tion and Bethe-Heitler pair production), synchrotron 
radiation and adiabatic expansion. The mesons pro- 
duced in photonmcson interactions always decay in 
astrophysical environments. However, they may suf- 
fer synchrotron losses before the decay, which is 
taken into account in this model. 

If the relativistic electrons are accelerated together 
with the protons at the same site, their injection 
spectrum shows most likely the same spectral shape 
cx j^''' with qe = qp. The relativistic primary e~ ra- 
diate synchrotron photons which constitute the low- 
energy bump in the blazar SED, and serve as the tar- 
get radiation field for proton-photon interactions and 
the pair-synchrotron cascade which subsequently de- 
velops. The SPB-model is designed for objects with 
a negligible external target photon component, and 
hence suitable for BL Lac objects. The cascade re- 
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Figure 5. Best-fit time-averaged leptonic (solid 
curves) and hadronic (dot-dashed curves) models for 
BL Lacertae in 2000 \Bdttcher & ReimeA \200A) . 



distributes the photon power to lower energies where 
the photons eventuaUy escape from the emission re- 
gion. The cascades can be initiated by photons from 
7r°-decay ("7r° cascade"), electrons from the tt"^ — * 
/i^ — > decay ("tt^ cascade"), p-synchrotron pho- 
tons ( "p-synchrotron cascade"), charged /z-, tt- and 
if-synchrotron photons ("/i='=-synchrotron cascade") 
and from proton-photon Bethe-Heitler pair pro- 
duction ("Bethe-Heitler cascade"). 

Because "tt" cascades" and "tt^ cascades" 
generate rather featureles s photon spectr a 
l|Mucke fc Protheroel l2001t iMucke et all I2003D . 
proton and muon synchrotron radiation and their 
reprocessed radiation turn out to be mainly re- 
sponsible for the high energy photon output in 
blazars. The contribution from the Bethe-Heitler 
cascades is mostly negligible. The low energy 
component is dominanted by synchrotron radiation 
from the primary e^, with a small contribution 
of synchrotron radiation from secondary electrons 
(produced by the p- and /^^-synchrotron cascade). 
A detailed description of the model itself, and its 
implementation as a (time-in dependent) Monte- 
Carlo code, h as been given in iMiicke fc Protherol 
l|200H) and .Reimer et alJ l|2004(l . This code has been 
used, e.g., to gener a te the spectral fits presented in 
iBottcher fc Reimed l|2004j) for BL Lacertae in 2000 
(see also Fig. ISJ. 



4. MODELING RESULTS FOR BL LACERTAE 
IN 2000 



Both the leptonic and hadronic models described 
in the previous sections have been applied success- 
fully to the data obtained during the multiwave- 
length observing carnpaign on BL Lacertae in 2000 
iBottchcr ct al., 200%. The best time- averaged 
spectral fits are shown in Figs. ^andlSl The best-fit 



Figure 6. Time- dependent model spectra from the 
best-fit leptonic model fo r BL Lacertae in 2000 
\Bdttcher & ReimeA. \200^) . The time sequence is: 
thin solid — > thin dotted — > thin long- dashed thin 
dot-dashed — > thin short-dashed — > thick solid — + 
thick dotted — > thick long-dashed, in equi-distant time 
steps of Aiobs ~ 



parameters for the leptonic and hadronic model fits 
differ substantially in the following ways: (a) The 
overall jet power required in the hadronic models is 
^ 2 orders of magnitude higher than for the leptonic 
models (~ 6 x 10'*'' ergs s~^ vs. <6 x lO''^ ergs s^^); 
(b) the magnetic fields in hadronic jet models are 
a factor of ~ 20 higher than for leptonic models 
(B - 30 - 40 G vs. 2 G), and the bulk Lorentz 
factor of the emission region is a factor of ^ 2 lower 
for the hadronic models (~ 7 - 9 vs. ~ 18). 

Considering the time-averaged emission of BL Lac- 
ertae in 2000, hadronic models predict a sustained 
level of multi-GeV - TeV emission which should 
be detectable with second-generation atmospheric 
Cherenkov telescope systems like VERITAS, HESS, 
or MAGIC. In contrast, our leptonic model only pre- 
dicts a peak fiux exceeding the anticipated nominal 
MAGIC sensitivity during short fiares; the accumu- 
lated fiuence over observing time scales of several 
hours might not be sufficient for a significant detec- 
tion. Thus, a future VHE detection of BL Lacertae 
would be a strong indication for hadronic processes 
being at work in this object. 

The detailed spectral evolution of our best-fit lep- 
tonic jet simulation is shown in Fig. |B| In order 
to determine the best-fit variability scenario, various 
generic fiaring scenarios had been investigated and 
compared to the observed optical and X-ray spectral 
variability patterns (see Fig. Specifically, we had 
investigated scenarios invoking a temporary increase 
in jet power, a flattening of the electron injection 
spectral index, an increasing high-energy cutoff 72 of 
the electron injection spectrum, and various combi- 
nation of these. Other parameter fluctuations could 
be ruled out on the basis of general, analytical con- 
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siderations, without detailed simulations. We found 
that the observed optical and X-ray spectral vari- 
ability in BL Lacertae in 2000 can be reproduced 
through short-term fluctuations of only the electron 
injection spectral index, with all other parameters 
remaining unchanged. 

Remarkably, our time-dependent fits indicated that 
an injection index larger than q 2.3, even dur- 
ing the peak of an individual short-term flare, is 
required. If the injection of ultrarelativistic elec- 
trons into the emitting volume is caused by Fermi 
acceleration at relativistic shocks, detailed numeri- 
cal studies have shown that with fully developed tur- 
bulence in the downstream region, a unique asymp- 
totic index of q ^ 2.2 2.3 should be expected 
(e.g., IXchterberg et all 12 001: Gallant ct al., 1999). 
However. recentlv lOstrowski fc Bednara 1,2002) have 
shown that Fermi acceleration might lead to drasti- 
cally steeper injection spectra if the turbulence is not 
fully developed. Furthermore, depending on the ori- 
entation of the magnetic field at the shock front, an 
abrupt steepening of the injection spectra may result 
if the shock transits from a subluminal to a super- 
luminal configuration. In this context, our leptonic 
fit results may indicate that such predominantly ge- 
ometric effects, may be the cause of the rapid vari- 
ability observed in BL Lacertae. 

Our spectral-variability simulation predicted 
counter-clockwise spectral hysteresis at X-ray ener- 
gies. Such hysteresis was not pre dicted in the specific 
SPB model fits presented in iBottcher fc Reirnerl 
l)2004j) . but could not clearly be ruled out either. 
Thus, sensitive spectral-hysteresis measurements 
of BL Lacertae could possibly serve as a test of 
our modeling results and a secondary diagnostic to 
distinguish between leptonic and hadronic models, 
though, by itself, it would not be sufficient as a 
model discriminant. 

iRavasio et al.l l)2003j) had previously noted the dis- 
crepancy bet ween the time-averag ed optical and X- 
ray spectra ijRavasio et all l2003(l which could not 
be joined smoothly by an absorbed power-law spec- 
trum. They had considered several possibilities to 
explain this discrepancy, including additional parti- 
cle populations, extreme Klein-Nishina effects on the 
electron cooling rates, and/or anomalies in the inter- 
galactic absorption. Our successful modeling of the 
observed time-dependent flux and hardness values at 
optical and X-ray frequencies in the framework of a 
leptonic model effectively removes the need for such 
additional assumptions and indicates that the time- 
averaging involved in compiling the detailed broad- 
band spectral energy distribution may be the cause 
of this apparent discrepancy. 
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